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Abstract. We consider second order elliptic divergence form systems with com- 
plex measurable coefficients A that are independent of the transversal coordinate, 
and prove that the set of A for which the boundary value problem with L2 Dirich- 
let or Neumann data is well posed, is an open set. Furthermore we prove that 
these boundary value problems are well posed when A is either Hermitean, block 
or constant. Our methods apply to more general systems of PDEs and as an ex- 
ample we prove perturbation results for boundary value problems for differential 
forms. 

MSG classes: 35J25, 35J55, 47N20 

1. Introduction 

We first review the situation for scalar equations. Consider the divergence form 
second order elliptic equation 

n 

(1) diYt,.A{x)Vt,,U{t, x) = J2 d^A^A^)d,U{t, x) = 0, 

i,j=0 

on the upper half space R^"^" := G R x R" ; t > 0}, n > 1, where the 

matrix A = (^ij(a^))"j=o ^ -^oo(R-"; -^(C^"^")) is assumed to be t-independent and 
strictly accretive with complex coefficients. In this generality, when no regularity is 
assumed of the coefficients, the natural conditions to impose on U at the boundary 
R" are one of the following. 

• Dirichlet problem (Dir-A): U{0,x) = u{x) for a given function u{x). 

• Neumann problem (Neu-A): Aoj{x)djU{0,x) = 4>{x), where is 
given. 

• Dirichlet regularity problem (Reg-A): diU{0,x) = diu{x), 1 < i < n, where 
u{x) is given. 

In this paper, we consider these boundary value problems (BVPs) in L2(R"), i.e. the 
boundary data are u G L2(R'^), (p G L2(R") and u G H^{TV^) respectively, and 
for well posedness a unique function U (t, x) with certain L2 estimates is required. 
Detailed definitions are given in Section O 

These BVPs arise naturally when considering BVPs for the Laplace equation 
on a Lipschitz domain Q in R". As the main problem here is a local one, the 
result for such domains can be derived from the scale invariant case of a Lipschitz 
graph domain, i.e. we assume that Q = {{t,x) ; t > g{x)} is the domain above 
the graph of some Lipschitz function g. Through a change of variables U{t,x) := 
V{t + g{x),x), an harmonic function V in Vl corresponds to U in R""*"^ satisfying 
([T]) with coefficients A = [1 -|- | Vxfi'P, — VxS'*; — VxS', /], and the respective boundary 
conditions carry over from dQ to R". The coeffi cents appearing from this pullback 
technique are referred to as being of Jacobian type, and are in particular real and 
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symmetric, as well as independent of the transversal coordinate t. In this case, 
solvability of the Dirichlet problem was first proved by Dahlberg ^TSj , and solvability 
of the Neumann and regularity problems was first proved by Jerison and Kenig [23] • 
Later Verchota [34] showed that these BVPs are solvable with the layer potential 
integral equation method. For general real symmetric matrices A, not being of the 
Jacobian type, well posedness of the Dirichlet problem was first proved by Jerison 
and Kenig ^24] , and the Neumann and regularity problems by Kenig and Pipher [2^ . 

It is natural to ask whether the BVPs for the fundamental elliptic equation ([1]) 
are well posed for more general coefficients. Obvious generalizations are coefficients 
A{t, x) with t-dependence, as well as more general non-symmetric or complex co- 
efficient matrices. In both cases, it is known that well posedness does not hold in 
general. Caffarelli, Fabes and Kenig [TT] observed that some regularity in the t- 
coordinate is necessary for well posedness, and Kenig, Koch, Pipher and Toro [26] 
and Kenig and Rule [28] gave examples where well posedness fail in any Lp for suf- 
ficiently non-symmetric, but t-independent real coefficients in the plane which are 
discontinuous at x = 0. However, on the positive side they show that, for given real 
non-symmetric coefficients in the plane, the Dirichlet problem is well posed in Lp for 
sufficiently large p, whereas the Neumann and regularity problems are well posed 
for p sufficiently close to 1. 

In this paper, we consider only t-independent coefficients, but allow on the other 
hand arbitrary complex, strictly accretive coefficients A G Loo(R"; /^(C^'*'")). As 
remarked, well posedness does not hold in general for the BVPs in this case. But 
our main result Theorem 12.21 shows that the sets of well posedness 



(2) WPiX) := {A ; (X-A) is well posed in L2(R")} C Loo(R"; i:(C^+")), 



where X denotes one of the three BVPs Dir, Neu or Reg, are all open sets. As 
discussed above, the sets of well posedness contain all real symmetric coefficients. 
Our theorem thus in particular shows well posedness for small complex perturbations 
of real symmetric coefficients. This has also been proved in [2J by Alfonseca, Auscher, 
Axelsson, Hofmann and Kim, with other methods using layer potential operators, 
and in [3] by Auscher, Axelsson and Hofmann. In fact, our methods here give 
the new result that well posedness holds for complex Hermitean matrices and their 
perturbations. 

One may ask what is the motivation for considering complex coefficients. However 
interesting it may be to know well posedness for complex matrices, a main motivation 
is that this feeds back to give perturbation estimates for real matrices. In fact, to 
show that the solution U varies continuously, for fixed boundary data, as A{x) 
varies continuously in Lqo within the subspace of real symmetric matrices, there is 
no known method which does not use bounds for complex BVPs. The observation 
being used is that bounds for complex BVPs imply analytic dependence on A and 
in particular Lipschitz regularity with respect to A G L^. 



Turning to other consequences of Theorem 12.21 well posedness is well known to 
hold for all constant coefficients A{x) = A, and our theorem thus yields well posed- 
ness for perturbations here as well. The Dirichlet problem was first shown to be well 
posed for small perturbations of constant matrices, by Fabes, Jerison and Kenig [20] , 
using the method of multilinear expansions. The Neumann and regularity problems 
are tackled in [2] and [3]. 
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It is also known that (Neu-A) and {Keg-A) are well posed in L2 for complex 
matrices of block form, i.e. such that Ag^i = = Ai o for all 1 < i < n. This is a 
non-trivial result and is in fact equivalent to the Kato square root problem, proved 
by Auscher, Hofmann, Lacey, Mcintosh and Tchamitchian [5]. Our theorem thus 
yields well posedness for small perturbations of complex block form matrices, and 
is in this sense a generalization of the Kato square root estimate. With the further 
assumption of pointwise resolvent kernel bounds, this result is also implicit in [2]. 
However, our methods in this paper require no such pointwise estimates. 

Let us now discuss the methods underlying Theorem 12. 2[ For the proof we use, 
following [3j, boundary equation methods involving a Cauchy operator E^. The 
name Cauchy operator is used since coincides with the Cauchy singular integral 
operator when A = I and n = 1. The first step of the proof is to rewrite the second 
order equation ([T]) as the equivalent first order system 

Uwt,^A{x)F{t,x) = 0, 
\^cnT\t^:rF{t,x) = 0, 

taking the gradient vector field F{t,x) := Vt,xU{t,x) as the unknown function in- 
stead of the potential U. The Cauchy operator E^ is related to in the same 
way that the classical Cauchy integral operator is related to the Cauchy-Riemann 
equations. Just as the Cauchy singular integral operator is a Fourier multiplier with 
symbol sgn(^), that is, belongs to the functional calculus of d/dx, the operator Ea 
belongs to the functional calculus of a first order differential operator T4. The bisec- 
torial operator —Ta in L2(R'^; C^"*"") is the infinitesimal generator for the system ([3]) 
in the sense that these equations are equivalent to dtF + T^F = 0. The fundamental 
problem is to prove that this operator T4 has a bounded holomorphic functional 
calculus, and as a consequence that the Cauchy operator Ea is bounded. Given 
this, the perturbation results for BVPs follow as a consequence. 

In [3], it was proved that \\Ea\\ < 00 when ||y4 — AqUoo < e and Aq is either 
real symmetric, block or constant. This paper made use of a rather lengthy per- 
turbation argument, and also used square function estimates of Dahlberg, Jerison 
and Kenig [16] and estimates of harmonic measure of Jerison and Kenig [21], for 
solutions to ([T]) in the real symmetric case. 

In this paper we prove the boundedness of the holomorphic functional calculus 
of Ta, for all complex A, directly from the quadratic estimates proved by Axelsson, 
Keith and Mcintosh |lOj. In this way, our results build on the proof of the Kato 
square root problem [3]. That \\Ea\\ < 00 for all complex A may come as a surprise, 
in view of the above mentioned counter-examples to well posedness of the BVPs 
for non-symmetric coefficients. However, it is important to note that Ea itself 
has nothing to do with BVPs, it is an infinitesimal generator associated with the 
differential equation, and is not related to the boundary conditions (except in the 
case of block form matrices). As a consequence of the boundedness of Ea = sgn(T4), 
we prove in Theorem 12.31 that there is a Hardy type splitting 

L2(R"; C^^") 3 f = -^'''Ir" + F Iru 

of boundary functions / into traces of satisfying ([3]) in R±^", with estimates 
II /II 2 ~ ||-P''''|r"||2 + II -^^ I R" II 2- That a BVP is well posed is the question whether 
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the full traces -F^Ir" of solutions to the equations in R""*"^ are in one-to-one correpon- 
dence with the normal components (Neumann problem) or tangential parts (Dirich- 
let regularity problem) respectively. It is this one-to-one correspondence which may 
fail for some complex A. What we prove here and use for the proof of Theorem 12.21 
is that the Hardy subspaces {-F"^|r"}, and the projections / t-^ -F^Ir" onto them, 
depend analytically on A. 

Finally we note that the methods developed here go beyond scalar elliptic equa- 
tions like ([T]). The natural framework here is rather EVPs for elliptic systems of 
partial differential equations, as should be clear from ([3]). Thus we shall formulate 
our results for divergence form elliptic systems of m second order equations (but 
the reader interested in scalar equations only can set m = 1 throughout). In this 
setting, our well posedness results are mostly new. Previously known results are 
limited to systems with coefficients of Jacobian type, or more generally constant 
coefficient systems on Lipschitz domains. 

Well posedness of (Dir-A), (Reg- A) and (Neu-A) with L2 boundary values have 
been obtained for the Stokes' system by Fabes, Kenig and Verchota [21], and of 
(Dir-y4) and (Reg- A) for the Lame system by Dahlberg, Kenig and Verchota [T8] . 
For general constant coefficient symmetric second order systems, solvability result 
for (Neu-A) and (Reg- A) are found in Fabes [I9]. Under the weaker Legendre- 
Hadamard ellipticity condition, (Dir-A) and (Reg- A) where solved by Gao ^22j. As 
for non-symmetic systems, Verchota and Vogel [35] obtained Lp solvability results for 
(DiT-A), (Reg- A) and (Neu-A) in the spirit of [2S] and [2H] for certain non-symmetric 
constant coefficient Legendre-Hadamard systems of two equations on polygons 
in the plane. For general elliptic systems, the Kato problem was solved by Auscher, 
Hofmann, Mcintosh and Tchamitchian |6j, a consequence being the well posedness 
of (Neu-y4) and (Reg-A) for elliptic systems with block form coefficients. 

Note that the puUback technique, from the Lipschitz domain fl to R""^^, described 
above, works for more general divergence form equations or systems. In this case, 
coefficients A{x) in Q are transformed into coefficients 



A{x) :-- 



1 -{VMx)y 
/ 



A{x) 



1 

-V^g{x) I 



in R^^^. We also remark that our methods are by no means limited to divergence 
form equations. In Section [HI we give solvability results for exterior differential 
systems (1221) for differential forms, as an example of this. This generalizes the first 
order system ([3]), which is the special case of ( l22l) for 1-forms. Furthermore, we note 
that time-harmonic Maxwell's equations on a Lipschitz domain can be written as a 
system of equations ( l22l) for 1 and 2-forms with lower order terms added, through 
the above puUback technique. Thus, although not directly applicable, (122!) is closely 
related to Maxwell's equations. Solvability of Maxwell's equations on Lipschitz 
domains is due to Mitrea [22], and more general EVPs for Dirac equations were 
solved by Mcintosh and Mitrea [3D]. In fact, the Cauchy integral boundary equation 
method used in this paper, as well as in [3], was developed for solving EVPs for 
Maxwell's and Dirac's equations in the PhD thesis [8] of the second author. Further 
elaborations of the ideas presented in this paper, along the lines of thought in [3] 
and [9], working with general inhomogeneous differential forms taking values in the 
full exterior algebra and allowing lower order terms, one should be able to extend 
the theory to cover both Dirac's and Maxwell's equations. 
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2. Notation and results 

Wc begin by giving the precise definition of well posedness of the EVPs discussed 
in the introduction, or rather the corresponding EVPs for systems. Throughout this 
paper, we use the notation X ^ Y and X <Y lor estimates to mean that there 
exists a constant C > 0, independent of the variables in the estimate, such that 
X/C < F < CX and X < CF, respectively. 

We write {eg, ci, . . . , e„} for the standard basis for R-'^+"' with cq "upward" pointing 
into R+'''", and write t = Xq for the vertical coordinate. For the vertical derivative, 
we write do — dt- For vectors v — (I'Dolf-oT' ^® write vq and f y for the normal 
and tangential parts of v, i.e. (vo)o = '^^o ('^o)f = when I < i < n, whereas 
(fii)" = vf when 1 < i < n and (f||)o = 0. Frequently, we shall identify normal 
vector fields v = vq with the corresponding scalar functions Vq = (^0)™=!- 

For systems, gradient and divergence act as {'Vt,xU)f = diU^ and {d\Vt,xF)^ ~ 
X^iLo^i-^i"' '^i^^ correponding tangential versions V^U — {yt,xU)\\ and (divj;^)" = 
Z]"=i diF^. With cm\t-xF = we understand that djF^ = diF°', for all i, j ~ 
0,...,n,a = l,...,m. Similarly, write curlj.F|| = if djF" = diF^, for all i, 
j = 1, . . . ,n,a = 1, . . . ,m. 

We consider divergence form second order elliptic systems 

n m 

(4) E E diA'lfi^)djU^it, ^) = 0, a = 1, . . . , m, 

i,j=Q 13=1 

on the half spaces R^"*"" := {(t,x) G R x R** ; ±i > 0}, n > 1, where the matrix 
A = {Aff{x))lfsQ-;:;^ e Loo(R";i:(C(i+")'")) is assumed to be t-independent with 
complex coefficients and strictly accretive on N (curly), in the sense that there exists 
K, > such that 

n m „ n m „ 

for all / e N(curl||) := {g G L2(R"; C(i+")") ; curl^(^||) = 0}. 
Equivalently, this means that a Carding inequality 

/ Re [A(diva;ii|| + VxUq), {diVxU^^ + 'VxUo))dx > k I (|divj;ii|||^ + \ VxUo\^)dx 
Jr" ^ ^ Jr" 
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holds for u = [mq, Mh] G ij^(R"; C'^^"'""')™), since diVa- has dense range in L2(R"; C™) 
and has dense range in {g E L2(R"; C""*) ; curl^.5f = 0}. Sphtting Qi'^+^)"^ into 
normal parts and tangential parts C"™", we write 



A{x)v 









f II 



It is then clear that implies that Aqq is pointwise strictly accretive and that Ay 
satisfies a strict Garding inequality 



(6) 



G R" 



The condition ([S]) lies between pointwise strict accretivity, i.e. 

(7) Re{A{x)v, v) > k\v\'^, for all v G C(^+")'" and a.e. x G R" 

and the R"^^ Garding inequality 



(8) 



/ / Re{A{x)Vt,x9{t, x),Vt,x9{t, x))dtdx > k, \Vt^g{t,x)\^dtdx 



for g G H^(R^^] C™'). Clearly, (JTj) implies ([5]), which in turn implies ([H]), as is seen 
by taking f{x) = Vt,xg{t,x) for fixed t and then integrating over t. Furthermore 
([8]), implies IQ, which is seen by taking g{t,x) := 'ip{et)uo{x) and integrating away 
t, for some ip G C^(R+). Letting e ^ oo and e — >• respectively proves ([U]). In fact, 
only the Garding inequality ioi g E Hq{II^^; C"^) is needed for this argument. 

When n = 1, ([5]) is equivalent to strong accretivity (171) since in this case has 
dense range in L2(R; C™) and N(curl||) = L2(R; C^™). On the other hand, if A is of 
block form, i.e. Aqu = A^^o = 0, then ([5]) is equivalent to (E]) and to the R""*"^ Garding 
inequality ([HD, for ifi(R!^+^ C™) as well as for H^iRl^^; C'"), since © implies ®. 
It is also known that the R"^^ Garding inequality (IHl) implies strong accretivity ([7]) 
when m = 1, so for scalar equations ([5]), ([8]) and ([7]) are all equivalent. 

On the functions U = (f/")^^^ satisfying (j3]), we impose one of the following 
boundary conditions. 

• (Dir-y4): f/°(0,x) = u°'{x) for a given function m G L2(R'^; C"). 

• (Neu-A): - 'Zj,(3Ao;j{x)djU'^{0,x) = where G L2(R"; C") is given. 

• (Reg- A): diU°{0,x) = diU°{x), l<i<n, where u E H\W; C") is given. 

The boundary value problems (Neu-A) and (Reg-A) can be viewed as problems 
concerning a first order partial differential system, and this is the point of view we 
take here. Indeed, consider the gradient vector fields 

F{t, x) = Vt,xU{t, x) : R^+" — > c'^^+n)m_ 

Since the scalar potentials U are in one-to-one correspondence with the curl-free 
vector fields F{t, x), modulo constants, we can take F rather than U as the unknown, 
and equation (j3]) for U is rewritten as the equivalent first order system for F. 
Since the coefficients A{x) are independent of t, it is natural to view F from the 
semigroup point of view F{t,x) = Ft{x) = Vt,xU{t,x) E C\R+; L2(R"; C(i+'^)'^)). 
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Definition 2.1. (i) We say that the boundary value problem (Neu-A) is well 
posed if for each boundary data G L2(R"; C™), there exists a unique func- 
tion 

F{t,x) = Ft{x) = Vt,,U{t,x) e Ci(R+;L2(R";C(i+")'")) 

which satisfies (E]) for t > 0, and has limits limj^oo Ft = and limj^o Ft = f 
in L2 norm, where the full boundary trace / satisfies the boundary condition 
— (A/)o = 0. More precisely, by F satisfying ([S]), we mean that dt^AF)^ = 
— div2.(AF)||, = VxFq and curl;j.F|| = 0, where dt is taken in the strong 
sense, and x-derivatives in the distributional sense. 

(ii) We say that the boundary value problem (Reg- A) is well posed if for each 
boundary data V^m G L2(R"'; C"'"), there exists a unique function F G 
C^(R+; L2(R"; 0(1+")™)) which satisfies © for t > 0, and has limits lim^^oo Ft = 
and limj^o Ft = f in L2 norm, where the full boundary trace / satisfies 
the boundary condition /y = VxU. 

(iii) The Dirichlet problem (Dir-A) is said to be well posed if for each u E 
L2(R"; C"^), there is a unique function 

Utix) = U{t,x) e Ci(R+;L2(R";C'")) 

such that VxU E CO(R+; L2(R"; C"'")), where U satisfies (il) for t > 0, 
limt^o Ut = u, limi^oo Ut = 0, limf^oo "^t^xUt = in L2 norm, and J^^ VxUs ds 
converges in L2 when ^ and ti — > 00. More precisely, by U satis- 
fying (HD, we mean that J^°°{{AVs,xUs)ii,^xv)ds = -{{AVt,xUt)o,v) for all 
V e Co°°(R";C"^). 

The Dirichlet problem (Dir-A) will also be rewritten as a BVP for the first order 
system However, here it is not appropriate to consider the gradient vector 

field Vt,xU, since the boundary condition is a condition on the potential U itself. 
Instead we use the point of view of harmonic conjugate functions, and write F = 
Ueo+F\\, where the tangential vector fields are conjugate functions in a generalized 
sense and F satisfies ([3]), which is viewed as a generalized Cauchy-Riemann system. 
Details of this are given in Lemma 1121 where it is shown that the Dirichlet problem 
(Dir-y4) for U is equivalent to an auxiliary Neumann problem (Neu"'"-y4) for F. 

Our main result, which we prove this in Section |U is the following. 

Theorem 2.2. The sets WP{Reg), WP{Neu) and WP{Dir), as defined in are 
all open subsets o/Loo(R"; >C(C*-"'^+"^'")). Each of the sets of well posedness contains 

(i) all Hermitean matrices A{x) = A{x)* (and in particular all real symmetric 
matrices), 

(ii) all block matrices where Ag'f (x) = = A1'q{x), 1 < i < n,l < a, (3 < m, 
and 

(iii) all constant matrices A{x) = A. 

The notion of well posedness used here departs from the standard variational 
one. However, we show in Section [5] that the solutions obtained here coincide with 
the solutions obtained through the Lax-Milgram Theorem when A belongs to the 
connected component of WP which contains /. This connected component includes 
the three classes (i), (ii) and (iii) specified in Theorem 12.21 The notion of well 
posedness used here coincides with that in [3j for the EVPs (Neu-A) and (Reg- A) . 
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However, for (Dir-A) the meaning of well posedness differs slightly from that in [3] , 
as we impose an extra integrability condition on V^?7 here. 

A natural function space for solutions Ft{x) e Ci(R+; L2(R"; C(^+")™)) to the 
EVPs is Loo(R+; L2), with norm sup^^g 11-^* lb- Two other norms which are important 
are the square function norm := \\Ft\\2t~^dt and the norm ||A^*(F)||2, 

using the modified non-tangential maximal function 

iV,(F)(x) :=supt-(i+^)/lF|U,(Q(i,.)), 

where Q{t, x) := [(1 — Co)t, (1 + Co)t] x B{x; Cit), for some fixed constants Cq G (0, 1), 
ci > 0. 

The key result underlying Theorem 12.21 which we prove in Section [3l is the fol- 
lowing result on Hardy type splittings of L2(R"'). 

Theorem 2.3. Let A E Loo(R"'; jC(C^"^+"'''")) be a t- independent, complex matrix 
function which is strictly accretive on A/(curl||). 

Then each f G A/(curl||) is in one-to-one correspondence with a pair of vector fields 
F^{x) = F^{t,x) G C\R±; L2{R!'; C(^+")™)) m R^+" satisfying ^ and having L2 
limits limf^o± Ft^ = f^ o-i^d lim^^-i-oo = 0, such that 

f = f^ + r- 

This splitting is topological, i.e. \\f\\2 ~ ||/^||2 + ||/~||2, and the vector fields 
satisfy norm equivalences 

||/||2~sup||Fi||2^ llltS^Filll ^ l|iV,(F)||2. 

t>0 

Moreover, the Hardy projections : / i— >• F^ = F^ depend locally Lipschitz con- 
tinuously on A in the sense that 

lli^i - FfjU < CP2 - A|U^(K.)||/||2, 

where C = C{kai, i^A2A\^i\\ooA\^2\\oo) o,nd where \\F\\x denotes any of the four 
norms above. 

Restricting our attention to the Dirichlet problem in the upper half space, Theo- 
rem 12.31 shows in particular the following. 

Corollary 2.4. Let A G Loo(R"; /^(C*^^"'""'^™)) be a t -independent, complex matrix 
function which is strictly accretive on A/( curly ) and assume that A G WP{Dir) . Then 
any function Ut{x) = U{t,x) G C^(R+; L2(R"; C™')) solving Q), with properties as 
in Definition \2. 1\ has estimates 

I \u\'^dx^ sup \Ut\'^dx^ / \N^{U)\^dx^ // \Vt.:,U\h dtdx, 

JR" t>0 jRn J Jin J Jr^ + '' 

where u = f/|R,n. If furthermore A is real (not necessarily symmetric) and m = 1, 
then Moser's local boundedness estimate [33] gives the pointwise estimate N^{U){x) ~ 
N^{U){x), where the standard non-tangential maximal function is N^{U){x) := 
^^V\y-x\<ct \Uit,y)\, for fixed < c < 00. 

Theorem 12.21 shows in particular that A G WP{Dir) if A is real symmetric. Even 
for real symmetric A our methods yield a new proof of the estimate between the 
square function and the non-tangential maximal function above, first proved by 
Dahlberg, Jerison and Kenig ^16j, using estimates of harmonic measure by Jerison 
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and Kenig [21], in the scalar case m = 1. In [3], these estimates were used to 
prove \\Ea\\ < oo for real symmetric A. Here we reverse the argument: we prove 
\\Ea\\ < oo independently and deduce from this the equivalences of norms. 

In the case m > 1 of systems, Dahlberg, Kenig, Pipher and Verchota [IT] have 
shown equivalence in Lp norm of the non-tangential maximal function and the square 
function, for general constant coefficient real symmetric systems on Lipschitz do- 
mains. Also, such equivalence has been shown for certain non-symmetric systems of 
two equations in the plane by Verchota and Vogel [36] . 



3. Cauchy operators and Hardy spaces 



The boundary equation methods for solving EVPs which are used here and in 
[3], are based on Cauchy operators Ea_, with associated Hardy type subspaces. In 
this section we prove quadratic estimates for Ea and deduce from this Theorem 12.31 
How these Cauchy operators are used to prove Theorem 12.21 is shown in Section H] 
and in particular Lemma 14.31 Note that the operators E^ themselves depend only 
on the differential system ([3]), and have nothing to do with the boundary conditions. 
Thus they are the same for both Neumann and Dirichlet problems. 

We start by rewriting the Equations ([3]) in terms of an "infinitesimal generator" 
Ta- Write v G 0^-^+"^"^ as f = [fcfy]*, where fo G C" and fy G C"*^, and introduce 
the auxiliary matrices 



A 



Aoo ^oi 
/ 



A :-- 



A 



110 





Aim 



if A 



Aoo 



A 

A, 



Oil 



in the normal/tangential splitting of C^^^"-*"*. Recall that A being strictly accretive 
on N (curly ), as in implies the accretivity estimates for the diagonal blocks Aqq 
and y4||||. Since Aqq is pointwise strictly accretive, it is invertible, and consequently 
A is invertible. This is not necessarily true for A. 

Splitting normal and tangential derivatives in ( 
equations is equivalent to 



we see that this system of 



dt{AF)o + dw,{AF), 
dtF. - V.Fo = 0, 



0, 



together with the constraint curl^-Fn = 0. Since 

(9) AF={AF)o + F^^ and AF = Fo + {AF)^^ 

we have shown that ([3]) is equivalent to 

'dt(AF)o + diYMF)ii = 0, 
dt(AF),-VMF)o = 0, 



together with the tangential constraint cur^Fn 
we get 



0. Combining the two equations. 



(10) 



dtF + TAF = 0, 



where Ta is the following operator. 



10 



PASCAL AUSCHER, ANDREAS AXELSSON, AND ALAN MCINTOSH 



div^. 
-V. 



m 



Definition 3.1. Let D be the self-adjoint differential operator D : 

L2(R"; C(i+-)-) with domain D{D) := {[/o, /„]* G L^iK^; C(i+«)-) ; V,/o, divj, G 
L2}. Define the infinitesimal generator for ([3]) to be the operator 

Aoo^(y4o||Vx + diVa;v4||o) AqqCHv^A 

-V,. 
with domain D{Ta) '■= Ar^D(D). Let the transformed coefficient matrix be 

-1 



Ta := A ^DA 



A := AA 



so that 

(11) Ta = a-\da)a. 

The following is the main algebraic result of the paper. Recall that N(curl||) = 
{g E L2(R"; C(i+")") ; curU(^||) = 0} and note that N(curl||) = R{D). 

Proposition 3.2. The transformed coefficient matrix A is bounded and strictly ac- 
cretive on A/(curl||), i.e. satisfies if and only if A is bounded and strictly accretive 

on A/(curl||). Moreover A = A. 

Proof. Assume that A is bounded and strictly accretive on N (curly). As noted above, 
A is invertible and thus A is bounded. Since A acts as identity on tangential vector 
fields, it is clear that A : N (curly) N (curly) is an isomorphism. Strict accretivity 
of A on N (curly) now follows from the formula 

Re{A(Af),Af) = Re{Af,Af) = Re ((/o, Aoo/o + AoJ,) + (Ayo/o + A,J„f,)) 
= Re ((Aoo/o + Aoii/ii, /o) + (^iio/o + = M^f, /)• 

The identity A = A is straightforward to verify, and this shows that the above 
argument is reversible. □ 

We are now in a position to analyze the operator Ta- Due to ( |TT|) . it suffices to 
study operators of the form DB in L2(R"'; C^), where D is a self-adjoint homoge- 
neous first order differential operator with constant coefficients, and _B is a bounded 
multiplication operator which is strictly accretive on R(-D), i.e. there exists k, > 
such that 

(12) Re{BDu, Du) > k\\Du\\2, for all u G D{D). 

The applications we have in mind are the specific operators D and B = A from 
Definition 13.11 in which case R(-D) = N (curly), as well as generalizations of these in 
Section [HI 

Define closed and open sectors and double sectors in the complex plane by 
:= {z eC ; I arg2;| < a;} U {0}, := U (-5*^+), 

5°+ ■.= {zeC ; z^O,\aTgz\< u}, S: := 5°+ U (-S,%), 
and define the angle of accretivity of B to be 

u := sup I arg(5/, /)| < n/2. 

/^0,/GR(D) 
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Proposition 3.3. Let D be a self-adjoint operator and let B be a bounded operator 
in L2(R"; C^) which satisfies Ulj. 

(i) The operator DB is a closed and densely defined uo-bisectorial operator, 
i.e. a{DB) C S^^, where u is the angle of accretivity of B. Moreover, there 
are resolvent bounds ||(A/ — DB)^^\\ < 1/dist (A, S^j) when X ^ S^^. 

(ii) The operator DB has range R{DB) = R{D) and null space N{DB) = 
B~^N{D), where B~^ denotes the inverse image, such that 

L2(R"; C^) = R{DB) © N{DB) 

topologically (but in general non-orthogonally) . 
(ill) The restriction of DB to R{D) = R{DB) is a closed and injective operator 
with dense range in R{D), with estimates on spectrum and resolvents as in 
(t). 

Proof. As a consequence of fll2p . it is verified tliat DB and B*D are closed and 
desely defined adjoint unbounded operators, and tlie topological splitting 

L2(R"; C™) = N(D) © B*R(D) 

follows from perturbing the orthogonal splitting L2 = N(D) © R{D) with B*, which 
satisfies (HSl) as well. Since R{DB) = M{B*D)^ = N{D)^ = R{D) and N{DB) = 
R{B*D)-^ = {B*R{D))-^, the splitting (ii) for DB follows, by stability of splittings 
under taking orthogonal complements. 

Since DB = {DB)\-^^(BO in the splitting (ii), it suffices to prove resolvent bounds 

for DB on R{D). To this end, let u G R{D) n D{DB) and / = (A/ - DB)u. Then 

Im {Bu, Am - DBu) = Im {Bu, /). 

Since D is self-adjoint, {Bu, DBu) G R and we get an estimate |Im (A(i?M, m))| < 
11/11, from which the resolvent bound < ||/||/dist (A, S*^) follows. □ 

These properties of closed operators of the form DB have been known for some 
time, see for example [IJ and [14J, at least in the case when B is strongly accretive. 
The following theorem has also been known for some time in the case when D is 
injective, as it derives from the special case D = —id/dx developed in [3T| (see also 
Lecture 8 of [1]). In this case DB is similar to the operator of differentiation on a 
Lipschitz graph 7, and the boundedness ofsgn{DB) is equivalent to the boundedness 
of the Cauchy singular integral on 7, proved originally by Calderon when B — I is 
sufficiently small [T2], and in general by Coifman, Mcintosh and Meyer [TSj . 

The proof of the following theorem however is more involved when D is not injec- 
tive. In the general case it was proved in flUl Theorem 3.1(iii)], building on results 
for the Kato problem by Auscher, Hofmann, Lacey, Mcintosh and Tchamitchian [5] . 
It is also possible to give a direct proof, as shown in [1]. 

Theorem 3.4. Let D be a self-adjoint homogeneous first order differential operator 
with constant coefficients such that 

P/II>I|V/|| forallfeR(D)r^D{D), 

and let B be a bounded multiplication operator in L2(R"; C^) which satisfies 



12 



PASCAL AUSCHER, ANDREAS AXELSSON, AND ALAN MCINTOSH 



(i) The operator DB satisfies quadratic estimates 

r°° fif 

/ \\tDB{l + [tDBfr'ft- ^ ll/f , for all f G R{D). 
Jo 1^ 

(ii) The operator DB has a bounded holomorphic functional calculus in R{D), 
i.e. for each hounded holomorphic function b{z) on a double sector S°, u < 
V < tt/2, the operator b{DB) in R{D) is bounded with estimates 

\\h{DB)U,^L,<\\h\\L^m- 

For the precise definition of operators b{DB) in the functional calculus of DB we 
refer to [r|. Note that the map Hoo{S°) 3 h ^ h{DB) G C{R{D)) is a continuous 
algebra homomorphism. 

We now return to the operator T4 of Definition I3.1[ Note that the isomorphism 
A in ( fTTl) maps the subspace R(-D) onto itself. 

Definition 3.5. Let Ti denote the closed subspace 

n := R(D) = N(curl||) = {/ G L2(R"; C(i+")"^) ; curl,/,, = 0}. 

of L2(R"; C(^+")'") 

In this notation, Proposition 13.31 and Theorem 13.41 with B = A, thus have the 
following corollary. 

Corollary 3.6. The operator Ta from Definition \3.1\ is an uj-bisectorial operator in 
L2(R"; C(i+")™), where uj is the angle of accretivity of A. Furthermore we have a 
splitting 

L2(R";C(^+")™) = R(7^©A/(T^)=7^©{[0,/|,]* ; div^A,,,,/,, = 0} 

in which Ta = T^ln © 0. The restriction of T^ to Ti is an injective operator with 
dense range in Ti, which satisfies quadratic estimates and has a bounded holomorphic 
functional calculus. 

Note that the restriction of to 7i, which we continue to denote by Ta, coincides 
with the operator Ta used in [3] for m = 1. 

Of importance to us are the following operators, which are bounded operators in 
7i because they belong to the functional calculus of Ta- 

• The characteristic functions 

±, . _ r 1 if ±Rez > 
X - I g if ± Re z < 

which give the generalised Hardy projections := x^(Ta). 

• The signum function sgn(2;) = x'^i^) ~ X (-2) which gives the generalised 
Cauchy operator Ea '■= sgn(T4) = E^ — E^- 

• The exponential functions e~*'^', t > 0, which give the operators e"*'"^-*'. Note 
that \z\ := zsgn{z) does not denote absolute value for non real z, but z h-^ \z\ 
is holomorphic on S°^2- 

Note that the quadratic estimates in Corollary 13.61 can be written as 

(13) mTA)frj^\\f\\\ feR{D), 
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where ip{z) := z/{l + z"^). The estimate < remains vahd for any holomorphic ip{z) 
on S° such that \ip{z)\ < min(|2;|", |^|~°) for some a > 0. If furthermore '^p\s;_^_ 7^ 
and il>\s°_ 7^ 0, then the estimate ~ holds. See p. 

Proof of Theorem \2.3[ By Corollary I3.6[ the infinitesimal generator has a bounded 
holomorphic functional calculus. With the notation introduced above, define Hardy 
type subspaces E^H := {E^f ; / G H}, so that 

7i = E^Ti © Ej^Ti. 

In terms of the operator T^, the vector fields / and F^ are related as 

f = f^ + f-^Eif = f^ = F±|r„ ^ e^*l^-l/± = F±. 

Indeed, differentiating F^ = e^^^'^^^f'^ at {t,x) E R^"", we have 

dtF^ = t\Ta\F^ = -T^e^*l^-^I(±EA/^) = -T^e^*l^-^l/± = -TaF^ 

since G E^Ti. Thus satisfies (JTOl) . which we have seen at the beginning of this 
section is equivalent to Conversely, each vector field F G C^(Il+] L2(R"'; C'-"'^"''"^''")) 
which satisfies ( II U I) for t > and has limits limt_^oo Ft = and lim^^o Ft = f in L2 
norm, is of the form F{t, x) = e~*'"^^'/ and / G E\l-L. To see this, split Ft = F^ + F^ 
where F^ G E^H. Applying E^ to (dUD gives the equations dtF^ ± IT^IF^^ = 0. 
Multiplying with suitable exponentials shows that, for fixed t > 0, e*-'^~*^'^'^'F/ is 
constant for s G (0,t) and e^^~^'^^'^^^F~ is constant for s G (t, cxd). Thus F^ = 
and Ft^ = e"*'"^-*'/, where / = lim^^^o-^t- The corresponding result for F~ is proved 
similarly. 

The equivalences of norms follow from the quadratic estimates for T4 as follows. 
Boundedness of the complementary projections E^ shows that ||/|| ~ ||/^|| + ||/~||, 
the uniform boundedness of e"*'"^-*', t > 0, shows that ||/^|| ~ sup_|_j>o ll-^i^^ll; and the 
quadratic estimates for shows that ^ \\f^\\- Finally ^ ||iV*(F=^)|| 

was proved in [31 Proposition 2.56] for m = 1. The extension to divergence form 
systems is straightforward. 

To verify Lipschitz continuity, one shows that the operators 

(14) / ^ 6(T^)/ -.n^n, 

(15) / ^ ib{tTA)f)tyo : n ^ L2((a, b);n), 

(16) / ^ {^{tTA)f)t>o ■■ n L2(R+, dt/t; H) 

depend analytically on A, where b and ip are bounded holomorphic function on S*^, 
ijj decays at and 00 as in (fT3ll . and the interval (a, b) is finite. One proceeds similar 
to fiUl Theorem 6.4] and [3l Lemma 2.41], starting from the analyticity of resolvents 
A I—*- [XI — Ta)~^ , using the quadratic estimates from Corollary 13.61 for the operator 
Ta and the fact the uniform limits of analytic functions are analytic. 

Lipschitz continuity can now be deduced from analyticity. Given Ai and A2 G 
Loo (R"; 0(1+")™) which are strictly accretive on 7i, define 

A{0■.= A, + aA2-A^)/\\A2-A^\\^ 

so that ( I— > A{() is analytic in a neighbourhood Q of the interval [0, \\A2 — Ai||oo]- 
Consider the analytic function ( ^-^ ^(C) ^ F^ = b{tTA{Q))f, where b{z) := 
e-*l^lx^(2), which has bounds ||F,^||2 < II/II2 in VL. Thus WdF^/dCh < ||/||2 in 
Vt, from which Lipschitz continuity of ^ 1— >^ follows, uniformly for all t > 0. 
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Lipschitz continuity of A for the square function norm and the norm of the 

non-tangential maximal function are proved similarly, using analyticity of (1161) . with 
ip{z) = ze~^^^x^{z)i ( |T5l) . with h{z) = e~*'^'x^(z), respectively. For the non- 
tangential maximal function, we refer to the proof of [3l Theorem 1.1] for further 
details. □ 



4. DiRICHLET AND NEUMANN BOUNDARY VALUE PROBLEMS 

In this section, we return to the Dirichlet and Neumann EVP's and use The- 
orem 12.31 to prove Theorem 12.21 We start by translating (Dir-A) to an auxiliary 
Neumann problem (Neu-'"-74), which consists in finding U solving (jlj) with boundary 
condition 

• Neumann problem (Neu-'--74): —dtU{0,x) = (p{x), where ip G L2(R"; C™") is 
given. 

More precisely, we use the following first order formulation of well posedness. 

Definition 4.1. We say that the boundary value problem (Neu"'"-y4) is well posed if 
for each boundary data ip G L2(R"; C'"), there exists a unique vector field 

F{t,x) = Ft{x) = Vt,xU{t,x) G Ci(R+;L2(R";C(^+")'")) 

which satisfies (EI) for t > 0, and has limits limj^oo Ft = and limj^o Ft = f in L2 
norm, where the full boundary trace / satisfies the boundary condition — /o = ip. 

Lemma 4.2. Given u = ~(p G L2(R"; C™'), we have a one-to-one correspondence 

/oo 
\/s,xUis,x) ds 

between solutions U{t, x) to (Dir-A) and solutions F{t,x) to (Neu^-A). In particular 
WP(Dir)= WP(Neu^), where 

WP(Neu^) := {A ; (Neu^-A) is well posed} C L^iYC; c(i+")"^). 

Proof. Assume that F solves (Neu"'"-y4) with boundary condition /o = u, and let 
U := Fq. Then Ut G C\R+; L2) and V^U = dtF^^ G C%K+;L2). The limits 
lim^^o Ut = u, lim^^oo Ut = and limf^oo ^t,xUt = are direct consequences of the 
limits limi^o Ft = f, limj^oo Ft = limf^oo tdtFt = 0, whereas 

j\xUsds= r dsF^^ds = F^^{ti)-F,{to) ^-f^^, (to,ti) — ^ (0,cx)). 

J to J to 

The function U satisfies (jl]) since 

/oo roo 
iVxV,{AVs,xUs)n)ds= / ds{VxV,{AF,)^^)ds 

= -{Vxv, (AFt),) = {v,diYx{AFt),) = -{v,dt{AFt)o) = -{v, {AVt,xUt)o). 

Conversely, assume that U solves (Dir-A) with boundary condition U\jin = —ip, 
and let F{t,x) := — Vs,xU{s,x) ds. This gives a well defined function since 
F = U - limt.^oo It' ^xUsds, and F G C\R+;L2). Clearly, limt^oo Ft = and 
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linif^o -^i = ~^ ~ h, where h := 'VxUgds. The vector field F satisfies since 
curU ];°° V,U{s, x) ds = 0, dtF^^ = V.U = V.Fq and 

/oo 
(AV,,,t/,)|| ds) 

/oo 



t)o) 



for all V e Co°°(R";C'"), so that dt{AFt)o = -div^.(AFi)||. This completes the 
proof. □ 

As in the proof of Theorem 12.31 we denote the upper Hardy type subspace of H 

by 

E+n = {f ; F e C\K^+'';n) solves 1^ in R^+", lim WFt-fh = lim ||Fj2 = 0}. 

Proof of Theorem \2.iA By Theorem 12.31 . the spectral projection E\ = x+{Ta) onto 
this subspace is bounded and depends Lipschitz continuously on A. We now observe 
that {Reg- A), (Neu-A) and (Neu-'--y4) are well posed if and only if 



E+n ^{ge L2(R"; C""^) ; curl,^ = 0} 
E+n L2(R";C"^) 

^2 



f ^ {Af)o, 
f^fo 



are isomorphisms respectively. Since E'^ depends continuously on A by Theorem 
the following lemma shows that the sets of well posedness are open. 

Lemma 4.3. Let Pt he hounded projections in a Hilhert space Ti which depend con- 
tinuously on a parameter t G (—5, 6), and let S : H ^ )C he a hounded operator into 
a Hilhert space /C. If S : PqH IC is an isomorphism, then there exists < e < 5, 
such that S : P^TY ^ K, is an isomorphism when \t\ < e. 

Proof. Consider the family of operators PqH 3 f ^ SPtf G /C between fixed spaces. 
By assumption and continuous dependence, they are invertible when |t| is small. 
Thus it suffices to prove that Pt : P^fK PtH is invertible when \t\ is small. This 
holds since (/ - Pq{Po - Pt))"^^o, Po{I - Pt{Pt - Po))"^ : PtH ^ PqH are seen to 
be left and right inverses respectively. □ 



What remains to be proved is that the three maps are isomorphisms when A is 
either Hermitean, block or constant. In fact, it suffices to prove this for (Reg- A) 
and (Neu-A), due to the following result proved in [31 Proposition 2.52]. 

Proposition 4.4. The houndary value prohlem (New^-A) is well posed if and only 
if (Reg- A*) is well posed. 

That (Reg-y4) and (Neu-A) are well posed for Hermitean, block and constant 
coefficents, follows from P, Section 3] when m = 1. For Hermitean and block form 
coefficients, these proofs are readily adapted to systems, but to be self contained, 
we give simplified proofs below. 
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Define the operator : = 



-1 
/ 



, which reflects a vector in R*^. Write A^"*" 




/ 



for the tangential projection and = |(1 — A^) 



the normal projection, so that A^ = A^''^ — A^ . 



1 




for 



4.1. Hermitean matrices. Let / G E\l-L. This means that there is a vector field 
Ft in Rij"" such that dt.Ft = —TAFt, limt^^ Ft = and lim^^o Ft = f- Recall that 



H—rOO -'■ t 

A'^DA 



and note that DN + ND = 0. Furthermore, assuming that A* = A, it is seen from 
the definition of A that the Hermitean condition translates to (A)* = NAN. The 
Rellich type identity which is useful here is the following. 



{NAf, Af) 



{NATAFt, AFt) + {NAFt, ATAFt)dt 



dt{NAFt,AFt)dt 



{NAD AFt, AFt) + {NAFt, DAFt)dt = / {{ND + DN)AFt, AFt)dt = 
'0 Jo 

Thus {{Af)o, {Af)o) = {{Af),, {Af),), or in view of ©, 

(17) ((/)o,(A/)o) = ((A/),„ (/)„). 

Consider first the Neumann problem. From ffT7|) it follows that 

ll/f ^ Re{Af, f) = Re {{Af)o, /o) + {{Af)„ /„)) = 2 Re((A/)o, /o) < ||(A/)o|| ||/||. 

This shows that ||/|| < ||(A/)o|| holds for the Neumann map E^H 3 f ^ (^/)o! 
which implies that this map is injective with closed range. 

It remains for us to prove surjectivity of this map. Note that the above estimates 
also show that E\p-i 3 f ^ {^tf)o is injective with closed range when At : = 
(1 — t)I + tA, < t < 1. It follows from the proof of Lemma [4.31 and the method of 
continuity that all these maps have the same index. Since / = € WP, it follows 
that A e WP. 

Well posedness for (Reg-A) and (Neu"'"-yl) is proved in a similar way, by keeping 
the factor /y and /o respectively from ffT7|l . 

4.2. Block matrices. Note that in this case the Neumann problems (Neu-A) and 
(Neu-'"-y4) coincide, and that the operator T4 has the form 



Ta 





-V. 



Aqq diVjjAii 




Note that in this case, the accretivity condition ^ splits into the two indepen- 
dent assumptions 'Re{AQQU,u) > ||m||2 and Re(74|||| V^f , V^f) ^ ||Vsf||| for all u G 
L2(R"; C") and V xV G L2(R"; C"*"). Since the diagonal elements in Ta are zero, 
so are the diagonal elements of Ea since 

Ea = ^a{^a) = {^a) 











where L := — Aoodiv^Aim V^:. Another way to see this is from the calculation 
NEa = Nsgn{TA)N-^N = sgn{NTAN-^)N = sgn{-TA)N = -EaN. 
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To prove well posedness, we need to prove that A^"*" : E^H N^H and : 
E^Ti. — >■ N~T-C are isomorphisms. From EaN + NEa = 0, we obtain explicit inverses 
as 



2Et : N+n — > Etn and 2Et : — > EtH. 



For example, to see that {2E\g) = g when g G N^Ti, we calculate 

2N+E+g = iV+(/ + EA)g = g + N+E^g 

= + i(l + N)EAg = g + \{EAg - EANg) = g. 

In fact well posedness of the Neumann and regularity problems for block coeffi- 
cients is equivalent to the Kato square root estimate 

II VLm||2 ~ ||V:i.m||2, 

as was first observed by Kenig [251 Remark 2.5.6]. To see this, we deduce from 
Equation that / G E;^n, i.e. / = Ea/, if and only if /n = -V^L-^/^fo, or 
inversely /o = L~^/^AoodiVa;y4||||/|| and curl^./!, = 0. This can be used to construct 
f = fo + Ea/o = Ea/w + f\\ from either /o or /y. Note that the Kato estimate 
translates to ||/o|| ~ ||/||||- 

4.3. Constant matrices. If A is a constant matrix, we can make use of the Fourier 
transform. 

(i) First consider the simpler case when m = 1. In this case, the solutions to the 
BVPs can be explicitly computed on the Fourier transform side, since the problem 
reduces to an eigenvector calculation. At the frequency point ^ 7^ 0, the space H 
corresponds to the two dimensional space := {zeo + ; z,w E C} of vectors 

" te] 



with tangential part parallel to ^. The operator D corresponds to : = 
Compressing the constant matrix A to Ti.^, we define the 2x2 strictly accretive 



< 



matrix 



a b 




' ^00 


Ami ' 




c d 




.e*^iio 




, so 



d 



Computing eigenvalues and vectors shows that zcq + G x±{.D^A) if and only if 
(19) 



{az + hw) = ^i(c — b) ± i\J ad — \{h + c) 



w. 



Applying the similarities in (fTTj) . we characterize well posedness as follows. That 
(Neu-A), (Reg-y4) and (Neu-'--A) are well posed means that zcq + G x±{.E)^A) is 
determined by z, w and az + bw respectively. This is straightforward to verify using 

(USD. 

(ii) Next consider the case m > 1. In this case, we perform a Rellich type argument 
on the Fourier symbol, or rather we make a "reverse Rellich estimate" . 

The space Ti^ is now isomorphic to C^*" since z,w E C". In view of Proposi- 
tion it suffices to prove a-priori estimates ||/|| < H/JI and ||/|| < ||(A/)o|| for 

/ G x+(T^)7Y| uniformly for almost all ^ G R", where := A D^A. However, 
since x+i^T^) = X+{T^) for t > 0, it suffices to consider the unit sphere |^| = 1. By 
continuity and compactness, we need only to verify that no non zero vector / such 
that /|| = or (A/)o = can be in the Hardy space, i.e. be of the form / = -F(O), 
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where F : R_|_ — > Ti^ satisfies dtF = —T^F and limj^oo -F = 0. To prove this, we use 



the fact that D| = / to obtain 



{D^Af,Af) = - / d,{D^AF{t),AF{t)) dt 
Jo 

POO POO 

= 2 ReiAF{t),AF{t))dt = 2 / Re{AF{t), F{t)) dt. 
Jo Jo 

We now observe that the left hand side vanishes if /y = or {Af)o = 0, and from the 
right hand side we then see that F = identically, and therefore / = 0. The method 
of continuity, perturbing A to I now shows that the maps / i-^ /y and / t— {Af)o 
are surjective, and thus isomorphisms. 

We have now completed the proof of Theorem 12.21 □ 



Proof of Corollary \2.4\ As in Lemma (4.21 a function U solving (jlj), with properties 



as in Definition 12.11 is the normal part of a vector field F = U + F^^ solving ([3]), 
with properties as in Definition 12. 1[ Theorem 12.31 shows that F = e"*'^-*'/, where 
f = F\-£in G E^Ti and that estimates 



^snp \\Fth^\\N,{F)h^\\\tdtFS 
t>o 

hold. If A e WP{Dtr) = WP{Neu^), then ||/||2 ^ \\u\\2 and ||Ft||2 ^ WUth for all 
if: > 0, since Ft G E^H. For the square function norm we observe that dtF = Vt^xU, 
and for the non-tangential maximal function clearly ||A^*(F)||2 > ||A^=k([/)||2 holds. 
As \\u\\2 ~ sup^>o \\Ut\\2 for solutions to (jlj) has been shown, we have \\Ut\\2 ^ \\Us\\2 
when t > s. The reverse estimate ||A^*(F)||2 < ||A^*(f/)||2 now follows from 

l|iV.(t/)f >sup / / / \U{s,y)\'dsdydx 

t>0 jR"J\y-x\<citJ\s-t\<CQt 



sup-/ llt/sTc^S > sup ||f/(i+co)tf ^ 

t>oJ\s-t\<cot t>0 

This proves the corollary. □ 



5. Uniqueness of solutions 

In this section we compare the solutions to the EVP's (Neu-A), (Dir-A) and (Reg- 
A) in the sense of Definition 12. ![ with the standard solutions obtained from the 
Lax-Milgram Theorem. This uses the homogeneous Sobolev space H^^R}^"", C™), 
equipped with the norm := /j^i+n \^t,xU\'^, and the subspace of functions with 

vanishing trace. Continuing our first order approach to EVP's via Q, we make the 
following definition. 

Definition 5.1. Introduce spaces of vector fields 

L^(R^+"; C(^+")™) := {F G L2(R^+"; C(^+")™) ; curlj^i+nF = 0} and 
L^°(R^+"; C(^+")™) := {F G i:2(R++"; C(^+")™) ; curlRi+n(F,) = 0}. 

The condition curlKi+n(F2) = here means that the extension by zero F^, of F 
to R^^", is curl free, or formally: curlj^i+nF = and the boundary trace of F 

is normal to R". If F G L^(R++"; 0^^+")'"), it is seen that there there exists 
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U e L2°'^(R^,_"^"; C™"), unique up to constants among the distributions on R^"^", such 
that 'Wt,xU = F. Define Hilbert spaces 

H\R]+''] C") := {U G 4°"(R^+"; C"*) ; Vt,^U G L^(RV'"; 0^^+")'")}, 

ifo'(R++"; C™) := {U G 4°"(R^+"; C") ; Vt,,f/ G L^"(R++"; C(^+")'")}, 

with norms so that the correpondence U ^ F = Vt,xU is an isometry. 

It is straightforward to verify that a function U G if ^(R^^"; C™) belongs to the 
subspace ffQ(Rj',.^"; C") if and only if there exists a constant C such that U extended 
by C to Ri+" belongs to if i(Ri+"; C"). 

Functions U G if^(R^,j'"; C™) are well defined only up to constants, whereas for 
U G ifg(Ri'j_^"; C™), we will choose the constant so that U\-r." = 0. It is not true 
that ifg(Ri',J^"; C™) C L2(R+'^"; C"*), as a scaling argument readily shows. However, 
Poincare's inequality shows that 



\U{t,x)\'—^^< II \Vt,xU{t,x)\'dtdx, f/Gifo^(R^+";C 



Ri+" 1 + r + \x\^ J jr,i^+" 



1 + t2 + |a;|2 
If F = Vt,xU solves ([3]), then we formally have 



Ja{U, V):= II {AVt.xU, Vt,xV) dtdx = - I ((A/)o, v) dx, 

J JR^+" JR" 

where / = i^|Rn and f = 1^|r". As pointed out in Section [21 the standing assumption 
that A is strictly accretive on N (curly), i.e. implies that the Garding inequality 
m in R^' holds, i.e. \JAiU,U)\ > \\U\\l,. 

Note that V ^-^ f-^„{(j){x),v{x)) dx in Equation (PUI) below defines a bounded func- 
tional on ifi(R^+";C'") if = div^.?^;, where ^ G L2(R";C'"), w G L2(R";C"™), 
since / |(/^(0P Hiax(|^|~^, l)d^ < oo and the trace map V t-^ v maps if-'^(Ry""; C"^) 
if ^/^(R"; C™). Furthermore V i-^ Jyi(PtM, V^) defines a bounded functional on 
ifoHR++"; C™) if M G ffi(R";C'") since JJ \^e-'\^\u{0?dtd^ ^ hWmiR^) < o^, 
where Pt denotes the Poisson extension 

p^^^),= r((l + ^)/2) f tu{y)dy 



7^(l+n)/2 y^,^ {t^ + \x- y|2)(l+n)/2 " 

The Lax-Milgram Theorem proves the existence and uniqueness of the following if ^ 
solutions U. 

Definition 5.2. We say that (p is good boundary data for (Neu-A) if L2(R"'; C™) 3 
(f) = divj..^, where w G i72(R'^; C"™). If is good, we define the solution to the 
Neumann problem to be the unique function U G if ^(R^""; C'") such that 

(20) JAiU, V)= I (0(x), v{x)) dx, for all V G ff ^(R^+"; C™). 

Jr" 

We say that u is g'ooc? boundary data for (Dir-A), or equivalently that VxU is g'ooc? 
boundary data for (Reg- A), if m G if ^(R"; C™'). If u is good, we define the if ^ solution 
to the Dirichlet (regularity) problem to be the unique function U G H^(R}^^; C™) 
such that 

(21) Ja{U, V) = 0, for all V G ifo(R++"; C'"), 
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and U{t,x) - Ptu{x) G i/o(R++"; C'"). 

The goal in this section is to prove the following uniqueness result. 

Theorem 5.3. Assume that A belongs to the connected component of WP(Neu) / 
WP(Reg) / WP(Dir) that contains I. If the boundary data is good, then the solutions 
to (Neu-A) / (Reg- A) / (Dir-A) in the sense of Definition \2. 1\ coincide with the 
solutions. 

Remark 5.4. For general A in the set of well posedness, the solutions constructed 
in this paper using the boundary equation method do not necessarily coincide with 
the solutions. Examples of this were shown in [7J. Note that these examples 
combined with Theorem 15.31 proves the existence of many coefficients that do not 
have well posed EVP's (even in when n = m = 1 with real A), sufficiently many to 
disconnect these A with non-if^ solutions from the identity. 

The proof of Theorem 15.31 uses the following lemma with Aq = I. 

Lemma 5.5. Let Aq be a block matrix. Then there exists e > 0, such that if 
11^ — ^olloo < e O'^d the boundary data is good, then the solutions to the EVP's 
in the sense of Definition \2.l[ coincide with the solutions. 

Proof for (Neu-A). Let F = Vt,xU = e~*'^^'/ be the boundary equation solution to 
(Neu-y4) in 11^^ with data = div^w = — (A/)o. Using the isomorphism A from 
(ITTi) . we define the similar Hardy function / := Af G x+{dA)TC. 

Let N~ and A^"*" be the normal and tangential projections from Section HI The 
boundary condition on / can be written N" f = —0. We solve for / by making the 
ansatz / = 2x+{DA)h, where h G N^H (i.e. /ly = 0). This yields the equation 

-0 = N-f = 2N-x+{DA)h = (/ + N-sgn{DA))h 

for h, in the normal subspace N'Ti. We note the following properties of "the 
double layer type operator" Ka '■= N~sgn{DA)N~ . (See [H] for explanations of 
this terminology.) When A = Aq, then as in Section 1^12) it follows that = 0, 
since the diagonal entries of sgn{DA) are zero in the normal/tangential splitting 
of the space. Theorem 12.31 shows that Ka depends continuously on A. Moreover 
Ka{R{D)) C R{D) since 

KADg = N-sgn{DA)DN+g = Dsgn{AD)N-^ g = D sgn{AD) g) 

when g G D{D). Therefore, when ||y4 — v4o||oo is small, we can expand (/ + Ka)~^ 
in a Neumann series and deduce that h G R(-D) since —0 = —Dw G R(-D). Indeed 

J2i-KAf{-Dw) = ^( - J2^~N+sgn{AD)N+)''w^ =: Dwn, 

k=0 k=0 

where and Dw^ h converges in L2. Since D is closed, h G R(-D). This shows 
that / = 2x+{DA)h G R(/^i), and thus / G R(Ta). In particular / = |Ta|^/Vo for 
some /o G H. Quadratic estimates for the operator Ta now shows that 

// \VtxU\'^dtdx= / \\Ftfdt 

JJYt}+^ ' Jo 

= ^°'ll(t|T^I)^/V*l^^'/off ^ll/of <oo. 
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Thus U e ifi(R^+"; C™). To verify let V e ifi(R^+"; C'") and consider the 
function 



git):= [ iiAFt)o,Vt)dx, t > 0, 



where we view t i-^ {AFt)^ as a C°° curve in R(div2.; L2) and t 1^ as a continuous 
curve in H^/^ilV; C™). If Vt G C\R+; H^/^), then 

g\t) = [ ((-div.AFi, Vt) + {AFu dtVt)) dx = f {AFt, Vt,.V) dx. 

Hence g(T) — g{e) = ff^^_^^rp{AFt,'Vt^xy)dtdx. This also holds for general V G 
H^(T{}^^] C"), which can be shown by mollifying t Vt. Taking limits (e, T) 
(0,00) proves (EOj). □ 

Proof for (Reg-A). Similar to the proof for the Neumann problem, we consider the 
equation 

V^M = N+f = 2N+x+{DA)h = {I + N+sgn{DA))h 

for h G N~^TC, in the tangential subspace. We deduce that the trace / of the solution 
is in the range of Ta, and therefore U G if^(R^+"; C™) and (EI]) follows as in the 
proof for the Neumann problem. 

To prove that U - PtU E H^iR}^''; C""), it suffices to show that doHj = djHo on 
R^+", for j = 1, . . . , n, when H is Vt,x{U — Ptu) extended by zero. To this end, let 
$ G C^(R^'^") and consider the function 



gj(t) := / Hj{t,x)^{t,x) dx, t>0. 
Jr" 



Since H is curl-free on R+^", we have dogj = J^n Ho{—dj^) + Hjdo^. As we have 
L2{IV^) convergence 

Hj = djUt — djPtu — > djU — djU = 0, t — > 0, 

integration of dogj over t G (0, 00) shows that J J HjOq^ = J J H^dj^ for all $ G 
Co°°(Ri+"), i.e. doHj = djHo. □ 

Proof for (Dir-A). Let F = Vt,xU be the boundary equation solution to (Neu-*-- 
A) with data ip = —u G i7^(R";C'^), and recall that ?7 = Fq is the solution to 
(Dir-y4). Consider the boundary trace / = F|pt" G E^H. We now instead use 
the isomorphism A from ( fTTj) . and define the similar Hardy function / := Af G 
Rix+{AD)). 

The boundary condition on / can be written A^~/ = N~ f = u. We solve for / 
using the ansatz / = 2x+{AD)h, where h G N'Ti. This yields the equation 

u = N-f = 2N~x+{AD)h = {I + N-sgn{AD))h 

for h. We note that the double layer type operator Ka '■= N^sgn{AD)N^ maps 
Ka{D{D)) C D{D) since 

DKAg = N+Dsgn{AD)N-g = N+sgn{DA)DN- g = N+sgn{DA)N+{Dg) 

when g G 0{D). As above, since \\Ka\\ is small when \\A — Ao||oo is small, we can 
expand (/ + Ka)~^ in Neumann series and deduce that h G D{D) since u G D(D). 
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This shows that / = 2x+{AD)h G D{AD), and thus / G D(Ta). In particular 
/ G D(|Tyi|^/^). Quadratic estimates for the operator T4 now shows that 

/» /• /»oo poo 

// \Vt,.U\^dtdx= / ||M||2dt= / llTAFtll^di 

J J-R^+" Jo Jo 

Jo 

Thus f/ G ifi(R^+"; C"^), and (EOl) follows as in the proof for the Neumann problem. 
Finally, note that Vt,xU = —T^Ft = — e~*'"^^'(T4/). This shows L2 convergence 

— ^ -(T4/)|| = VJo = VxU. 

As in the proof for the regularity problem, U — PtU G Hq(TV^"; C") follows. □ 

Remark 5.6. (i) Note that for any L2 boundary data, the solution to the Neumann 
and regularity problem always satisfies WVt^xUW^ dt < 00, whereas the solution 
to the Dirichlet problem always satisfies WVt^xUW^ dt < 00. Thus the problem 
whether good boundary data give solutions concerns large t for the Neumann 
and regularity problem, and small t for the Dirichlet problem. 

(ii) The structure of the problem in Lemma 15.51 is best explained abstractly as 
follows. Let Til ^ 'Ho be a continuous and dense inclusion of Hilbert spaces. Assume 
that To : Tio Tio is an isomorphism which restricts to a bounded operator Ti : 
Hi — *• Til. It follows from [29i Theorem 11.1] that we have regularity Hi = Tq^(Hi) 
if and only if Ti is a Fredholm operator and has index zero. See also [H Proposition 
3.2.16]. 

In our situation, Tq = I + K, Hq = N^H and Hi is either R{D) n N^H or 
D{D) n N~H. In principle, the technique of Lemma 15.51 could be used to prove 
regularity for more general A in the component of WP containing J. The problem 
though is that in general the well posedness of two different EVP's, for the matrix 
A, is needed both for the proof that Tq is an isomorphism and that Ti is Fredholm. 
Index zero for Ti could then be proved by the method of continuity, perturbing A 
to /. 

Proof of Theorem \5.3[ Fix good boundary data and let < a < 6 < 00. For all A 
with the assumed properties, let F = Fa = Vt,xU be the solutions given by The- 
orem [2]2l and let F^ = F^ = Vt,xU denote the standard solutions constructed 
with the Lax-Milgram Theorem as in Definition 12. 1[ 

From the quadratic estimates for T4 it follows with arguments as in [3^, proof of 
Theorem 1.1] that A ^ F4 G L2(R" x (a, 6); C(^+")'") is analytic on WP. The main 
result this uses is the analyticity oi A^ b{TA) for operators &(Ta) in the functional 
calculus of Ta- This was proved in [TUl Theorem 6.4]. Moreover, it is straightforward 
to verify that A F^ E L2{'R^^^; is analytic. This means that whenever 

Aq G WP, C D D E z ^ A{z) are coefficients depending analytically on a complex 
variable z, and A(0) = Aq and H G i^2(R" x («, b); c(^+")'"), then the scalar function 
z I— >• h{A{z)), where 

h{A) := [[ (Fa - F^, H) dtdx, 

J J{a,b)x'R" 

is analytic on D. 



SOLVABILITY OF ELLIPTIC SYSTEMS 



23 



Consider one of the EVP's and fix A in the connected component of WP containing 
/. Pick a sequence of balls = B^Ak] r^) C WP, k = 0,1, . . . , N, such that Aq = I, 
A^- = A and S^-i n5fc 7^ 0. We may take tq < e, so that h = on Bq hj Lemma [^7S1 
Now assume that h = on Bk-i and pick any A^ G Bk- Let A^ G -Bfc-i fl Bk and 
let A{z) := (1 — z)A^ + zA^. Then h{A{z)) vanishes on a neighbourhood of 0. By 
analytic continuation h{A^) = h{A{l)) = 0, and since A^ G Bk was arbitrary, h = 
on Bk- We arrive at the conclusion that h{A) = 0. Since a,b and H are arbitrary, 
it follows that Fa = F^. □ 



6. Boundary value problems for differential forms 



In this section, we demonstrate how Theorem 12.21 and Theorem 12.31 generalize to 
exterior /interior differential systems for fc-vector fields, i.e. differential forms of order 
k. 

We use the notation from [3l Section 2.1]. In particular, for fixed A; G {1, 2, ... , n}, 
we consider functions 



F{t,x) = ^F{s,,...,sk}{t,x) Cs, 



A e 



taking values in the space = a'^R^^" of complex fc- vectors on R^+". The vector 
fields in ([3]) is the special case k = 1. We point out that we assume the component 
functions Fg to be scalar valued here (m=l), although the methods apply, mutatis 
mutandis, to systems of exterior differential systems. A natural generalization of 
the first order system is the interior /exterior differential system 



dl^BF{t, x) = 0, 
dt,^F{t,x)=0, 



(22) 

where F : R^^" — > A*'. Here the exterior and interior derivative operators are 

n 

dt,,F = Vt,, A F = ^ A djF = fidtF + d,F, 



j=0 



dLF 



t,x 



J2 J djF = -fx*dtF + dlF, 

j=0 



where a denotes exterior product and j denotes (left) interior product, and /i/ = 
Co A / and /i*/ = Cq j /. The matrix function B G Loo(R"; £(a'^)) is assumed to be 
t-independent and pointwise strictly accretive in the sense that 

Re{B{x)w,w) > /tlwp, for all w G a''' and a.e. x G R". 

To prove an analogue of Theorem 12.31 for the Equation ( l22i) . we proceed as in 
Section [3] and introduce auxiliary matrices 



B :-- 



B±± B^ 
/ 



B :-- 



B, 





5,M 



if 5 



5xx 

B,, I 



B^ 

Bn 



A e 



in the normal/ tangential splitting of A'^. Recall that a basis /c-vector a 
is normal if one of the factors is Cq, and tangential otherwise. Denote the tangential 
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and normal parts of / by /y and f^. Splitting each of the Equations (1221) into normal 
and tangential parts and using the analogue of ©, shows that (!22|) is equivalent to 

(dt(BF),-fidl{BF), = 0, 
\dt{BF), + fi*dMF')^ = 0, 

together with the constraints d^F^^ = = d*{BF)j^. These tangential derivatives in 
the equations define the appropriate function space 

■■= {/ e L2(R"; a'^) ; = = dUBf),} 

generalizing Ti from Definition 13.51 Note that when k > 2, the space depends 
on B, unlike the case k = 1. 

The normal derivatives in the equation give an equation dtF + TbF = 0, where 
the infinitesimal generator is 

(23) Tb := B'^DB. 

Here D := ^*dx — fid^ is a self-adjoint differential operator. The operator Tb has 
similarities 

b~\db)b = Tb = b-\bd)b, 

where B := B_B ^ is shown to be pointwise strictly accretive as in Proposition 13. 2[ 
Thus Proposition 13.31 applies and proves that Tb = B {DB)B is an cu-bisectorial 
operator, uo being the angle of accretivity of B. Moreover, Tb restricts to an injective 
tu-bisectorial operator in = B R(-D) with dense range, and in the splitting 

L2(R'^;A'=) = R(T^©N(Tb)= 7^1 ©{[/„/„]* ; 4/x = = 

we have Tb = Tgl^fc © 0. 

Similar to the proof of Theorem 12.31 Theorem 13.41 proves the boundedness of the 
Cauchy operator 

Eb ■■= sgn (Tb) 

and the Hardy projections E^ := x±(^b)- To handle perturbation theory for the 
variable space Hb, we extend this operator to L2(R"; A'^) by defining E^f = Eb/ = 
when / G N(Tb). 

We obtain the following result on Hardy space splittings of H% C L2(R"; A'^). 

Theorem 6.1. Let B G Loo(R"; >C(a'^)) be a t -independent, complex coefficient 
matrix function which is pointwise strictly accretive. 

Then each f G is in one-to-one correspondence with a pair of k-vector fields 
F^ = T±(t, ■) G Ci(R±;L2(R"; A'^)) m R^+" satisfying ^ and having L2 limits 
lim(_>o± = ^.''T'd lim^^-i-oo = 0, such that 

f = r + r- 

Under this correspondence, we have equivalences of norms II/II2 ~ ll/^lb + 11/ II2 
and 

(24) ||/±||2^sup||F±||2^|||ta*F±|||. 

±t>0 

Moreover, the Hardy space projections L2(R"; A'^) 9 / t-^ F^ = Eg := e^^^'^^^E^f 
depend locally Lipschitz continuously on B in the sense that 

- F|JU < C7||i?2 - i?l|U^(R.)||/||2, 
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where C = C{k,Bi, i^B2i ||-Bi||oo7 II-B2II00) o^^^^ where \\F\\x denotes any of the three 
norms in 



We remark that the proof of the non-tangential maximal estimate ||A^*(F)||2 ~ 
1 2 from [3j Proposition 2.56] in Theorem 12.31 uses the divergence form structure of 
the second order system. This technique does not generalize to more general exterior 
differential systems. 

Finally we extend the results in Section H] and show how Theorem 16.11 gives per- 
turbation results for EVP's for fc-vector fields. The natural EVP's are the following. 
We are looking for a A;-vector field Ft G Ci(R+; L2(R"; A*=)) solving ([22D in 
with L2 limits limf^o+ Ft = f and hmt^oo Ft = 0, where the boundary trace / 
satisfies one of the following. 

• Tangential boundary condition [Tan-B): /y = g, where the given boundary 
data g G L2(R"; A'^) is tangential and satisfies d^g = 0. 

• Conormal boundary condition (Nor-i?): {Bf)^ = g, where the given bound- 
ary data g G L2(R"; A'^) is normal and satisfies d^g = 0. 

Note that when k = 1, (Tan-i?) coincides with the Dirichlet regularity problem 
(Reg-i?) and (Nor-i?) coincides with the Neumann problem (Neu-i?). 

Theorem 6.2. The sets of well posedness WP{Tan) and WP{Nor) are both open 
subsets 0/ Loo(R"; i3(A^)) and each contains 

(i) all Hermitean matrices B{x) = B{x)* (and in particular all real symmetric 
matrices), 

(ii) all block matrices B(x) = ^^^^^ ^ 
im] all constant matrices Bix) = B. 



, and 



What is new here as compared with [3], is the perturbation result around Her- 
mitean and constant matrices, as well as the openness of the sets of well posedness. 
The proof of Theorem 16.21 is similar to the proof Theorem 12.21 We observe that 
(Tan-i?) and (Hoi-B) are well posed if and only if 

^B^l ^{g^ ^2(R"; a'^) ; ii*g = 0, 4^7 = 0} : / ^ /,„ 

^B^l -^{g^ ^2(R^ a'^) ■ ^^g = Q,dlg = Q}■.fv-^ {Bf),, 

are isomorphisms respectively. Theorem 16.11 shows that Eb depends continuously 
on B, so we obtain from Lemma [4.31 that WP(Tan) and WP(Nor) are open sets. 

That Hermitean and block matrices belong to WP is shown as in Sections 14.11 and 
14.21 mutatis mutandis. For constant matrices, the reverse Rellich argument used for 
second order divergence form elliptic systems m > 1 applies. For exterior differential 
systems, the symbol of the operator is 

:= B ^ D^B_, := i{lJ.*H( + 

acting in the 2 (^"^) -dimensional space := {/ G A'' ; fifi^f = = fi*fi'^Bf}. The 
proof uses that if fifi^f = = ii*^i*^f , i.e. / G R{D^), then 

(25) Oy = -{fi* fi^fifil + fifilfi*fi^)f = {fi^fi*fifil + filfin*fi^)f 
For the anticommutation relations, we refer to [H Lemma 2.3]. 
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